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ABSTRACT 

The discrepancy between cosmological Li abundance inferred from Population II dwarf 
stars and that derived from Big Bang nucleosynthesis calculations is still far from being 
satisfactorily solved. We investigated, as an alternative route, the use of Li abundances 
in Population II lower red giant branch stars as empirical diagnostic of the cosmological 
Li. Both theory and observations suggest that the surface Li abundance in metal poor 
red giants after the completion of the first dredge-up and before the red giant branch 
bump, arc significantly less sensitive to the efficiency of atomic diffusion, compared 
with dwarf stars. The surface Li abundances in these objects - after the dilution 
caused by the first dredge-up - are predicted to be sensitive to the total Li content 
left in the star, i.e. they are affected only by the total amount of Li eventually burned 
during the previous main sequence phase. Standard stellar models computed under 
different physical assumptions show that the inclusion of the atomic diffusion has an 
impact of about 0.07 dex in the determination of the primordial Li abundance - much 
smaller than the case of metal poor main sequence-turn off stars - and it is basically 
unaffected by reasonable variations other parameters (overshooting, age, initial He 
abundance, mixing length). We have determined from spectroscopy the surface Li 
content of 17 Halo lower red giant branch stars, in the metallicity range between 
[Fe/H]~ —3.4 and ~ —1.4 dex, evolving before the extra-mixing episode that sets 
in at the red giant branch bump. The initial Li (customarily taken as estimate of 
the cosmological Li abundance A(Li)o) has then been inferred by accounting for the 
difference between initial and post-dredge up Li abundances in the appropriate stellar 
models. It depends mainly on the T e // scale adopted in the spectroscopic analysis, 
and is only weakly sensitive to the efficiency of atomic diffusion in the models, so long 
as one neglects Li destruction caused by the process competing with atomic diffusion. 
Our final A(Li)o estimate spans a relatively narrow range, between 2.28 and 2.46 dex, 
and is ^0.3-0.4 dex lower than predictions from Big Bang nucleosynthesis calculations. 
These values of A(Li)o are corroborated by the analysis of samples of red giants in 
the Galactic globular clusters NGC 6397, NGC 6752 and M4. Our result provides an 
independent quantitative estimate of the difference with the Big Bang value, and sets 
a very robust constraint for the physical processes invoked to resolve this discrepancy. 
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ulation II - (Galaxy:) globular clusters: individual (M4, NGC6397, NGC6752) 



1 INTRODUCTION 



* Based on observations collected at the ESO-VLT under pro- 
grams 068.D-0546, 072.B-0585, 266.D-5655, 077.D-0018, 65.L- 
0165 and on data available in the ELODIE archive. These re- 
search has also made use of the SIMBAD database, operated at 
CDS, Strasbourg, France. 



The discovery by ISpite fc Spite] (|l982al lbT) of a uniform 
Li abundance in the atmospheres of main sequence (MS) 
Halo field stars with T eff above ~5500 - 5900 K and 
[Fe/H] below about —1.5 dex - the so-called Spite 
Plateau - has been widely interpreted as a signature 
of the cosmological abundance of lithium, produced dur- 
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ing the Big Bang nucleosynthesis (BBN). Adopting the 
standard notation A(Li)=log(n(Li)/n(H)) + 12, the Spite 
Plateau value is between 2.1 and 2.4, depending mostly 



on th e adopted T e ff scale ( s ee, e. g., Bonifacio fe Molard 



19971; ICharbonnel fc Primasl 120051; lAsplund et al 



2006 



Bonifacio et al.ll2007l ; I Aoki et al.ll2009l ; ISbordone et aljboiC , 

and references therein). 

Estimates of the cosmological baryon density obtained 
from the power spect rum of the cosmic microwave back- 
ground fluctuations (|Spergel et al.l I2007T ) combined with 
standard BBN calculat ions, predict howev er a cosmologi- 
cal A(Li) =2.72±0.06 (|Cvburt et alj 120081 ). Taken at face 
value, the BBN A(Li) is higher - by at least a factor of 2 - 
than A(Li) measured in Spite Plateau stars. 

The interpretation of Spite Plateau abundances in 
terms of the initial chemical composition of Halo stars needs 
however to take into account the possible effect of stellar 
evolution on the star surface chemical composition. Stellar 
models that include only convection as element transport 
mechanism predict for Spite Plateau stars a negligible pre- 
MS Li depletion, and no changes of the surface abundances 
during the whole MS phase. With these assumptions, the 
observed A(Li) is essentially equal to the primordial value. 

On the other hand, detailed stellar evolution calcula- 
tions need to include additional transport mechanisms that 
follow from first principles. In low mass stars, atomic diffu- 
sion - i.e. the slow transport of chemicals due to tempera- 
ture, pressure and abundance gradients - is the main process 
that needs to be included, for taking diffusion into account, 
the agreement between solar models and co nstraints from 
helio seismology is greatly improved (see, i.e.. iBahcall et all 
Il997h . A general prediction of models including atomic dif- 
fusion is that Population II MS stars in the Spite Plateau re- 
gion, with their shallow convective envelopes, are efficiently 
depleted of Li an d metals over the MS li fetime (see, e.g., the 
seminal paper bv lDelivannis et al.lll990h . Moreover, the sur- 
face Li depletion is a function of the star initial metallicity 
(and T e ff), because of the associated variation of the depth 
of the convective envelopes. These predictions are at odds 
with the observed constant (both as a functi on of [Fe/H] and 
T e f f) A(Li) along the plateau (but see also ISbordone et al.l 
l20ld . for the evidences of a meltdown of the Spite plateau in 
the very metal-poor r egime) . Inclusion of radiative levitation 
(|Richard et al.ll2005l ) leaves the overall picture unchanged. 
For [Fe/H] below ~ —2.0 dex radiative levitation moder- 
ates the surface depletion of Li due to atomic diffusion but 
the general picture is not qualitatively modified, whereas at 
higher metallicities the effect on the surface Li is negligible 
0. Overall, the depletion of Li in Spite Plateau stars due to 
atomic diffusion (plus radiative levitation) can reach several 
tenths of a dex at the lowest metallicity and/or higher T e //. 

For the models with diffusion (and radiative levita- 
tion) to display uniform A(Li) abundances along the Spite 
Plateau, one needs to include some ad-hoc turbulent mix- 
ing in the radiative zone bel ow the convec t ive b oundary, 
that limits the settling of Li. iRichard et al. (2005) include 
a specific parametrization of turbulent mixing, treated as a 
diffusive process, with a free parameter that fixes the values 



1 One has to notice that in I Richard et all I ll2005h the term atomic 
diffusion includes also the process of radiative levitation 



of the turbulent diffusion coefficient. A suitable calibration 
of this parameter - that enables also the transport of addi- 
tional lithium down to burning temperatures - produces a 
flat Li abundance along the plateau, with A(Li) depleted by 
a bout 0.4 dex with re spect to the initial value (see Fig. 8 
in lRichard et al.ll2005h . This depletion, when applied to the 
range of A(Li) measured by different authors, can poten- 
tially bring into agreement the Spite Plateau with BBN pre- 
dictions. 

IPiaul I 2008) includes instea d the effect of a rotation- 
ally induced tachocline mixing (|Spiegel fc Zahnl Il992l ) as 
employed in solar models (|Brun et al.lll999T ) to improve the 
agreement between theoretical and observed sound speed. 
Also in this case mixing is treated as a diffusive process, 
and the author discusses the calibration of the associated 
free parameters and the dependence on the assumed rotation 
history. As a result, tachocline mixing - that in plateau stars 
moderates the efficiency of atomic diffusion - appears to be 
able to produce a uniform A(Li) along the Spite Plateau, 
with a depletion of 0.2 dex compared to the initial value. 
Given the range of measured A(Li), the agreement with the 
BBN lithium abundance is still marginal at best. 

Similar estimates of primordial Li based on spec- 
troscopy of MS stars in globular clusters are complicated 
by the fact that in many cases the observed A(Li) display 
spreads r elated to the presenc e of second generatio n stars 
(see, e.g.. lBonifacio et al.ll2007l ; iD'Orazi et alj|20icl . for the 
case of 47 Tuc). The metal poor cluster NGC6397 is one 
object with appar ently an ess e ntially homogeneous initial 
Li abundance and iKorn et al.l ( 2006a ) spectr o scopic mea- 
surements plus the models by Richard et al.l (|2005l 'l with 
a suitably calibrated efficiency of their turbulent mixing 
parametrization - that turns out to be different from the 
case of field halo dwarfs - are able to reproduce the ob- 
served A(Li) (and [Fe/H]) measurements at the MS turn 
off and the base of the red giant branch (RGB), for an ini- 
tial A(Li)=2.54. The models fail however to match the ob- 
served post dredge-up A(Li) on the lower RGB, befo re the 
RGB bump. A recent reanalysis bv lLind et al.l (|2009l ) have 
revised down this value, to A(Li)=2.46. A detailed compar- 
ison with the observe d post dredge-up Li abundance is not 
presented. Moreover, iLind et "all |2009) conclude: "We find 
that some turbulence, in a very limited efficiency-range, is 
indeed required to explain observations. These models" (in- 
cludi ng atomic diffusion plus turbulent mixing parametrized 
as in lRichard et al.ll2005l ) "fail to reproduce the behaviour of 
Li abundance with effective temperature along the plateau, 
suggesting that a detailed understanding of the physics re- 
sponsible for depletion is still lacking." Note also that the 
assumption of a different temperature scale leads to different 
results, in terms o f match between observatio n s and models, 
as pointed out bv lGonzalez Hernandez et al. I (|2009D . 

Another cluster that does not sho w signatures of a 
spread in the initial Li abundance is M4. iMucciarelli et all 
(|201ll ) shows that models including diffusion plus turbulence 
with the same efficiency invoked for NGC6397, fail in case of 
M4. It appears that, in order to reproduce the Li abundances 
from the turn off to the lower RGB and satisfy at the same 
time the BBN constraint, for this cluster the turbulence has 
to reach deeper regions compared to NGC6397, region hot 
enough for some additional Li burning to occur. The trend 
of lithium abundance with effective temperature along the 
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sub-giant branch, found i n M4 by iMucciarelli et al. 
is similar to that found bv lGonzalez Hernandez et al 
in NGC6397. 

This brief summary of previous investigations shows 
how problematic, from a theoretical point of view, it is to de- 
termine the primordial Li abundance of Spite Plateau stars. 
In this paper we investigate a complementary avenue to es- 
timate the Halo primordial A (Li). We employ spectroscopy 
of Hal o stars evolving along the lower RGB, defined - fol- 
lowing iGratton et~alT l^OOof ) - as the portion of the RGB 
brighter than the luminosity corresponding to the comple- 
tion of the first dredge-up and fainter than the RGB-bump. 
We will show that the effect of atomic diffusion on the sur- 
face Li abundances of these objects is much smaller than for 
Spite Plateau stars. As a consequence, the observed A(Li) 
can be employed to set independent, very strong constraints 
on any additional physical process (i.e., turbulent mixing, 
preprocessing of Halo material, modifications to the BBN) 
eventually needed to reconcile these values with BBN pre- 
dictions. 

The paper is structured as follows. Section 2 analyzes 
the theoretical advantages of employing A(Li) measured in 
the atmosphere of lower RGB Halo stars to estimate their 
primordial Li abundance, while Sect. 3 discusses our selected 
sample of lower RGB halo objects and the derivation of 
A(Li) and [Fe/H]. Section 4 presents theoretical predictions 
for the surface Li abundance in Population II lower RGB 
stars - from models including atomic diffusion and convec- 
tion as element transport mechanisms - and the derivation 
of the initial Li for our selected star sample, that provides 
an estimate of the Halo primordial Li abundance indepen- 
dent of the Spite Plateau. Section 5 discusses the Li abun- 
dance in the lower RGB of the Galactic globulars NGC 6397, 
NGC 6752 and M4. A summary and conclusions close the 
paper. 



2 WHY LOWER RGB STARS ? 

The surface Li abundance after the first dredge-up is essen- 
tially a consequence of the dilution due to the increased size 
of the convective envelope after the MS turn off (plus a minor 
contribution from Li burning in the deep layers of the fully 
mixed envelope in the most metal rich Halo stars). At the 
end of the MS phase, when the deepening convective bound- 
ary reaches layers where the Li-burning (Tf,„ rn ~2.5 10 6 K) 
was efficient during the MS, the surface A(Li) begins to de- 
crease. This depletion essentially ends when the convective 
envelope attains its maximum depth and the first dredge-up 
is complete. An important point to notice is that atomic dif- 
fusion during the MS produces a local maximum of the Li 
abundance in the r adiative layers right below the convective 
envelope (see, e.g.. lRichard et al]|2005l ). for only a relatively 
small fraction of the envelope Li is transported deep enough 
to reach layers where it is eventually burned. As a conse- 
quence, models without and with diffusion - fully efficient 
or even moderated by levitation or some turbulence that 
just mixes back material diffused from the envelope - dilute 
similar amount of Li within the deepening convective region. 
In addition, the maximum size of the convective envelope is 
also weakly affected by diffusion - this can be also inferred 
by the fact that the predicted RGB bump luminosities with 



and without diffusio n are very similar |Cassisi et al.|[l997l ; 
iMichaud et al.|[2010T ) - and the resulting A(Li) abundances 
on the lower RGB are only slightly changed. This is in con- 
trast with models for upper MS stars, where the effect of dif- 
fusion on the surface abundan ces can reach several te nths of 
dex (see for example Fig. 3 in IMucciarelli et al.ll201lh . The- 
oretical models also predict that along the RGB, after the 
completion of the dredge-up, atomic diffusion and levitation 
are not able to modify app reciably the surface abunda nce of 
Li and all other elements jMichaud et al.ll2007i . |2010| ). 

The chemica l abun d ance measu r emen ts in giant stars 
bvlGratton et all (|2000r i, ISpite et all (120051 ) and lLind et~al] 
|2009l ). suggest that any additional element transport along 
the RGB is very likely inefficient in this phase. On the other 
hand, stars evolved beyond the RGB bump display the ef- 
fect of an additional mixing event, for which the so-called 
'thermohaline mixing' is nowadays the most popular candi- 
date (see, e.g. ICharbonnel fe Lagardd l2010l and references 
therein). 

These considerations lead to the conclusion that in gen- 
eral there is a simpler relationship between initial and cur- 
rent surface Li abundances in lower RGB stars compared 
to Spite Plateau objects. The basic reason is that, for a 
fixed initial Li abundance, lower RGB surface abundances 
are sensitive to the total amount of Li left in the star. In 
fact, after the effect of dilution is accounted for, they are 
affected only by the total amount of Li burned during the 
MS phase due to atomic diffusion plus possible additional 
element transport mechanisms - if they are invoked to solve 
the discrepancy with BBN results. On the other hand, the 
observed abundances along the Spite Plateau are determined 
by the evolution of the rate of Li depletion from the convec- 
tive envelope due just to diffusion or moderated/enhanced 
by additional turbulent processes. Combining the informa- 
tion obtained from Li in lower RGB stars with Spite Plateau 
data, sets much stronger constraints on how these proposed 
additional mixing processes act between the base of the con- 
vective envelope and the deeper Li burning regions. 

In this study, we address the issue of the Halo primor- 
dial Li abundance by coupling predictions from standard 
RGB stellar models to the Li abundances measured on a 
sample of lower RGB Population II field stars. As discussed 
before, this will provide an independent quantitative esti- 
mate of the difference (if any) with the Big Bang value, and 
set a very robust constraint on the efficiency of additional 
physical processes invoked to resolve this discrepancy, com- 
plementary to the interpretation of the abundances along 
the Spite Plateau. If the solution of the discrepancy involves 
pre-processing during the Galaxy formation or modification s 
to the BBN (see, e.g.. |Piau et~atll2006l ; ICvburt et alJl201Ch . 
our analysis will provide a solid quantitative estimate of, 
respectively, the net amount of Li burned during the early 
Galaxy evolution, or exactly how much Li must be synthe- 
sized in the 'revised' BBN. 

To the best of our knowledge, analyses of the Li abun- 
dances measured along the RGB have been so far mainly 
aimed at testing the agreement between observed and pre- 
dicted Li depletion after the first dredge-up, assuming a 
value for the initi al A(Li), often equal to wha t measured 
along the plateau (I Garcia Perez fc Primadljxjh . 
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3 LI ABUNDANCES IN LOWER RGB FIELD 
HALO STARS 

This section presents the observational data employed in our 
study, and the abundance analysis. Particular care is paid to 
the issue of the effective temperature scale and its influence 
on the derived values of A(Li). 



3.1 Dataset 

We have selected a sample of 17 metal-poor stars (with 
[Fe/H]<-1 dex), located in the part of the RGB fainter 
than the RGB bump magnitude level. The selection 
of this sample was done by cros s che cking data by 
ICavrel de Strobel, Soubiran fc Ralitel (|200lh . the SIMBAD 
website E|, previous works on chemical abundances in metal- 
poor giant st ars |Burris et al.l |2000| ; iGratton et al.l |2000| ; 
I Johnson! [20021) and the ESO and ELODIE archives, in or- 
der to select suitable high-resolution spectra, covering the 
Li line at 6708 A. 

Spectra of 11 stars w ere retrieved from th e ESO archive 
and observed with UVES (|Dekker et alJlioOtil ). employing si- 
multaneously different gratings. We used spectra obtained 
with the CD#3 cross-disperser with a nominal spectral res- 
olution of —54000 and a spectral coverage between —4760 
and 6820 A, to measure Li and Fe abundances, and spectra 
obtained with the grating CD#1 (-3050-3870 A) or CD#2 
(-3760-4980 A) to measure C abundances and the 12 C/ 13 C 
isotopic ratios. For 6 stars we u sed spectra from the ELODIE 
archive (|Moultaka et alj|2004h . with a spectral resolution of 
—42000 and a coverage between —4000 and 6800 A. 

Table [1] summarizes the main information for these 
stars (identification number of the corresponding catalog, 
adopted instrument, atmospherical parameters, abundances 
and SNR per pixel around the Li doublet). 



3.2 Abundance analysis 

The abundance analysis was performed by mean s of the 
codes SYNT HE (to compute synthetic spectra, IKurucj 
Il993bl . 120051 ) and WIDTH (that predicts the abundance of 
a given species b y matching the observed and theoretical 
equivalent widths, Kurucz 1993a; Castclli 2005), coupled to 
the ATLAS9 model atmospheres (lKuruczlll993al. |2005[). We 
used the Linux version of all the codes |Sbordone et al.ll2004l ; 
ISbordonj2005t ). The ATLAS9 models employed in this work 
were com puted with the new set of Opacity Distribution 
Function (|Castelli fe Kurucdl2003h and without the inclu- 
sion of the approximate overshooting in the calculation of 
the convective flux. 

Iron abundances were obtained from line equivalent 
widths (EWs) usi ng WIDTH; EWs were me asured using the 
DAOSPEC code (stetson fe Pancind |200S| ) . that performs 
an automatic measurement of the EWs under the gaussian 
profile approximation. The adopted linelist includes a hun- 
dredth of neutral iron lines whose oscillator strengths are 
from the recen t comp ilation of accurate laboratory log gf by 
iFuhr fe Wiesel (|2006f ) and —15-20 single ionized iron lines 



bv lRaassen fe Uvlingsl (|l99ST ) . For each star we refine itera- 
tively the linelist in order to include only lines predicted to 
be unblended from the inspection of synthetic spectra com- 
puted with the atmospherical parameters and metallicity of 
each target. The resulted iron abundances were s caled to the 
solar iron abundance of 7.52 |Caffau et al.ll2010T ). Total un- 
certainties in the Fe abundances were computed by taking 
into account the internal error (obtained as the dispersion 
of the mean normalized to the root mean square of the num- 
ber of used lines) and the errors related to the atmospherical 
parameters. 

Li abundances were derived with the help of spectral 
synthesis, in order to take into account the hyperfine struc- 
ture and the isotopic splitting of the Li resonance dou- 
blet at 6707.8 A. T he linelist for the Li transitions is from 
lYan fe Drake! (1 19951). Correcti ons for departures from LTE 
are from Carlsson et al.1 (| 19941 '). The total error in the A(Li) 
abundance is computed by taking into account the two main 
source of uncertainty: (i) the error in the adopted T e //, typ- 
ically 5A{Li)/8T eff -0.09-0.10 dex per 100 K for a giant 
star (giant stars are slightly more sensitive to T e // with 
respect to the dwarf stars, that are sensitive at a level of 
SA(Li)/ST eff -0.06-0.07 dex per 100 K); (ii) the error in 
the fitting procedure, obtained by employing MonteCarlo 
simulations: for each star, we injected Poisson noise into the 
best-fit synthetic spectrum (in order to reproduce the ob- 
served SNR around the Li line) and we repeated the analysis 
for 1000 MonteCarlo events. We estimated a la level from 
the resulting abundance distributions; due to the high SNR 
of the analyzed spectra (>100), the typical errors are lower 
than 0.03 dex. The errors due to uncertainties on gravity 
and microturbulent velocity are negligible (of the order of 
0.01 dex or less). 

We measured also the isotopic ratio 12 C/ 13 C in order to 
robustly assess whether the stars have experienced the ad- 
ditional mixing at/after the RGB bump, and establish their 
precise evolutionary stage. Being the spectra obtained with 
different instruments and configurations, we need to use dif- 
ferent indicators for the C abundance and its isotopic ratio. 
For the ELODIE spectra and for those stars observed with 
the UVES grating CD#2, C abundance and 12 C/ 13 C* ratio 
were estimated by spectral synthesis of the 2 A- 2 II band of 
CH (the G band) at -4310 A. For the stars for which the 
grating CD#2 is not available, we used the bluest portion of 
the spectrum (obtained with the CD#1 grism), covering the 
spectral region of 3000-3800 A. For these stars C was ob- 
tained from the 2 E- 2 n band of CH at 3143 A and from other 
CH transitions in the range between 3088 and 3115 A. These 
regions have been used also to measure the 12 C/ 13 C ratio. 
The atomic and molecular linelists were taken from the Ku- 
rucz compilation^]. For the CH transitions, the log gf by the 
Kurucz database were revised downward by 0.3 dex, in order 
to wel l reproduce the solar-flux s pectrum by Neckel fe Labsl 
|l984l ) with the C abundance bv lCaffau et al.1 ^Old ); simi- 
lar corrections were a pplie d also by other autho rs, see e.g.. 
iBonifacio et~al] (|l998r i and lLucatello et al.1 (|2003T ). As a first 
guess we derived C by assuming the solar 12 C/ 13 C ratio; 
then, we computed synthetic spectra varying the 12 C/ 13 C 
ratio, keeping C fixed. The procedure was repeated until 
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the convergence within a tolerance of 0.1 dex in the C abun- 
dance. Only upper limits can be inferred for several stars in 
our sample, due to the low intensity of the 13 CH transitions 
(only for the most metal-poor star in the sample, namely 
CD -30°0298, we cannot provide a reliable and strict upper 
limit). 

3.3 Lithium abundance from excitation 
temperatures 

We performed a fully spectroscopic analysis of the target 
stars, to derive their atmospherical parameters. They were 
computed by imposing no trend between neutral iron lines 
abundances and excitation potential x (to constrain T e ff), 
no trend between neutral iron lines abundances and the 
reduced equivalent width log EW/X (to constrain the mi- 
croturbulent velocity) and the same abundance (within the 
quoted uncertainties) from neutral and single ionized iron 
lines (to constrain the gravity). 

We noted that the iron lines with x <~1 eV exhibit 
abundances systematically higher than the other lines for 
stars with metallicity lower than -2.5 dex. Even if the 
adopted T e ff provides no dependence between Fe I abun- 
dance and x m the range of x ~1 _ 5 eV, the low-^ lines 
give abundances higher by ~0.3 dex. We excluded from the 
analysis the Fe I lines with x <1 e V, because their inclu- 
sion can force a too low T e // (with a reduction of about 
200 K), in order to erase a spurious slope. For more metal- 
rich stars, this discrepancy is reduced or totally erased. 
This effect was already noted in other studies of metal- 
poor stars (iNorris.Rvan fe Beersll200ll ; ICarretta et~ai]|2002l ; 
ICavrel et al. 2004 ) and probably ascribable to some inade- 
quacies of the model atmospheres based on the assumptions 
of 1-dimensional geometry and/or LTE: in fact, the low-x 
lines are typically formed in the outermost layers of the pho- 
tosphere, that are exposed to the UV radiation coming from 
the deep layers (thus, the non-LTE effects are more pro- 
nounced in metal-poor stars because of the low opacity and 
high transparency of the photosphere) . 

The derived iron abundances range from [Fe/H]= -3.40 
dex to [Fe/H]= -1.43 dex. The average Li abundance in- 
ferred from these 17 giants turns out to be A(Li)= 0.97 dex 
(a= 0.06 dex). 

Figure [1] shows the behaviour of the sample stars as a 
function of [Fe/H] and T e ff (upper and lower panel, respec- 
tively). We checked for the occurrence of slopes in these two 
planes, performing least-square fits by conside ring the er- 
rors i n both variables (following the approach bv lPress et alj 
Il992h . In the A(Li)-[Fe/H] plane we derived a slope of 
0.018 dex/dex (with an error of 0.044 dex/dex), while in 
the A(U)-T eff we find a slope of 0.013 dex/100 K (with an 
error of 0.028 dex/100 K). 

Eight stars in our sample are in common with 
iGratton etail (|2000h . We find a difference in the adopted 



temperatures of T ef f-T^ff 



-107 K (cr= 41 K), prob- 



ably due to their T e // scale, based on the grids of synthetic 
colors computed by R. L. Kurucz by means of the ATLAS9 
model atmospheres with the inclusion of the approximate 
overshooting and adopting the old set of Opacity Distri- 
butio n Function (we refer the reader to ICastelli fe Kurucz] 
120031 . for details about the differences between the two set of 
ATLAS9 models). The mean A(Li) difference for the stars in 
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Figure 1. Behaviour of the Li abundance as a function of the iron 
abundance (upper panel) and of the temperature (lower panel). 
The dashed line shown in the upper panel represents the average 
A(Li) of the field stars of the sample. 



common is A(Li)-A(Li) Gra " on00 = -0.19 dex (a= 0.07 dex), 
partially due to the different T e /f scales. The residual dis- 
crepancy most likely stems from the adoption of overshoot- 
ing model atmosp heres, that provide higher Li abund ances, 

as pointed out bv lMolaro. Bonifacio fe Primasl dl995h. 

T en stars are in common with Garcia Perez fe Primasl 
(|2006h : we find average differences T ef f-T^j cia - Pereis06 = - 
80 K ( ff = 84 K) and A(Li)-A(Li) Garcia - Perez06 = -0.11 dex 



(a= 0.11 dex). 

Figure [2] shows the position of the target stars in the 
T e ff -log g plane , toget her with theoretical isochrones by 
iPietrinferni et all (|2006h . computed with an age of 12.5 Gyr 
and metallicities Z=0.0001, 0.0003, 0.0006 and 0.001. These 
isochrones were computed with the same input physics, 
code, metal distribution and AY/AZ of our calculations. 
The position of the RGB bump is marked for reference as 
a grey shaded region. All the stars are located (within the 
uncertainties in the atmospherical parameters) below the 
RGB bump level, as confirmed also by the 12 C/ 13 C ratio 
higher than ~15-20. In fact, the occurrence of the extra- 
mixing episode after the RGB bump decreases dramatically 
the 12 C/ 13 C ratio, reaching values lower tha n 10 (see for in- 
stance \GtM^^^^^^M^^Q00^). This confirms, 
together with the homogeneous distribution of A(Li), that 
all these stars have not yet experienced the extra-mixing 
episode. The figure also shows that essentially all stars in 
our sample have attained the final surface A(Li) values, be- 
fore the onset of the post RGB bump extra mixing. 



3.4 Lithium abundance from Infrared Flux 
Method temperatures 

As a sanity check to assess the robustness of the average 
A(Li) value we repeated the analysis of the target stars 
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Table 1. Identification numbers, adopted instrument, temperatures, gravities, microturbu- 
lent velocities, [Fe/H] and A(Li) abundances, 12 C/ 13 C isotopic ratio and SNR around the 
Li line, for the observed stars. 



ID 


Instrument 


T„ ff 


lo£ £ 


Hi 


[Fe/Hl 

L / J 


AfLi) 


12 c/ 13 c 


SNR 






(K) 




(km/s) 


(dex) 


(dex) 






HD 2665 


T~1T /MM T — 1 

ELOD1E 


4950 


2.20 


1.80 


-2.18 


0.95 


10 


250 


HD 6755 


ELODIE 


5100 


2.30 


1.50 


-1.57 


1.01 


>20 


110 


HD 26169 


UVES 


5050 


1.90 


1.50 


-2.43 


0.99 


>30 


320 


HD 27928 


UVES 


5000 


2.00 


1.50 


-2.33 


1.04 


>30 


370 


HD 45282 


ELODIE 


5200 


2.80 


1.50 


-1.56 


1.01 


>30 


230 


HD 87140 


ELODIE 


5050 


2.30 


1.50 


-1.83 


0.92 


35 


150 


HD 111721 


UVES 


4950 


2.30 


1.30 


-1.43 


0.98 


40 


370 


HD 128279 


UVES 


5300 


2.60 


1.70 


-2.21 


1.02 


>35 


280 


HD 175305 


ELODIE 


5050 


2.40 


1.50 


-1.45 


1.07 


30 


180 


HD 200654 


UVES 


5150 


2.10 


1.50 


-2.96 


0.98 


>35 


290 


HD 211998 


UVES 


5200 


2.90 


1.20 


-1.59 


0.97 


40 


540 


HD 218857 


UVES 


5000 


1.90 


1.30 


-2.05 


0.92 


>25 


330 


HD 274939 


UVES 


5000 


2.00 


1.50 


-1.74 


0.96 


>20 


310 


BD -01°2582 


UVES 


5000 


1.90 


1.10 


-2.39 


0.82 


40 


280 


BD +23°3130 


ELODIE 


5200 


2.20 


1.50 


-2.59 


0.92 


>15 


170 


CD -24° 1782 


UVES 


5000 


2.10 


0.80 


-2.91 


0.96 


>20 


280 


CD -30° 0298 


UVES 


5100 


2.10 


1.80 


-3.40 


0.99 




330 



with other T e // scales, the temperature being the most 
crucial parameter in the derivation of A(Li). We infer T e // 
in our targets by means of suitable transformations between 
dereddened broad-band colors and effective temperatures 
obtained through the classical In f rared Flux Method 
(IRFM; iBlackwell. Petford fc Shallis! Il980l ). Several T e / / 
scales based on this technique are available in literature (e.g. 
Montegriffo et all 1 19981; [Al onso. Arribas fc Martinez-Rogerl 



Ramirez fc Melendej 120051; 

Gonzalez Hernandez fc Bonifacid 120091 ; ICasagrande et all 



19991; 



201(1 ). that make use of different model atmospheres, pho- 



tometric systems and recipes for the absolute calibrations 
and zero-points. The comparison between the various IRFM 
scales is beyond the purpose of this work and we focus our 
at tention only on the widely used in literat ure calibration 
bv lAlonso. Arribas fc Martinez-Rogerl ll 19991. hereafter A99) 
and t he recent one by iGonzalez Hernandez fc Bonifacid 
|2009l . hereafter G09), based on the Two Micron All Sky 
Survey 2MASS photometric data to compute the infrared 
monochromatic flux of the targets. 

In order to have a homogeneous set of magnitudes, 
we adopted for all the field stars the J and K s magni- 
tudes available in the final release of the 2MASS catalogue 
jSkrutskie et alj|2006l ). Note that the 2MASS magnitudes of 
the targets were transformed to the photometric system of 
the Telescope Carlos Sanchez (where the calibration b y A99 
is d efined) by means of the relations b y Carpente 
and lAlonso. Arribas fc Martinez-Rogerl (| 19981 ). 

The (J — A')o colors wer e obtai n ed by adopting 
the extinction coefficients by iMcCalll |200J) and the 
color excess E(B-V) from the i nfrared dust maps by 
ISchlegel. Finkbeiner fc Davis! (I1998T). corrected fol lowing the 



prescriptions by Bonifacio. Monai fc Beersi (|2000l ). 

For some stars, E(B-V) is larger than 0.4 mag and the 
derived T e f f is incompatible with the evolutionary stage of 
the targets and with the spectroscopic T e f f . For these stars 
we computed the color excess E(B-V) from the EW of the 



interstellar double t Na I lines (at 5890 and 5 896 A), through 
the calibration bv lMunari fc Zwitteil (| 19971 ). 

Gravities and microturbulent velocities have been de- 
rived spectroscopically, as described in the previous Section. 

The differences between the photometric and spectro- 
scopic scales are equal to T e ( jf-T^ c = +102 K (cr= 55 



k) and i, Ky ir'ff 



+4 K (<t= 57 K). This difference is 
fully consistent with the intrinsic difference between the two 
scales (see the discussion in G09). Consequently, we derived 
an average A(Li) of 0.97 dex (<r= 0.07 dex) and 1.07 dex 
(cr= 0.07 dex) when T e // by A99 and G09 are adopted, re- 
spectively. 

Thus, the spectroscopic T e // and those obtained with 
the A99 calibration can be considered on the same scale, 
while the scale by G09 is slightly hotter, providing Li abun- 
dances ~0.1 dex higher. 

Finally, we note that the general agreement between 
spectroscopic and photometric T e f j (apart from the different 
zero-points) seems to indicate that no relevant departures 
from the LTE condition occurs, at least when the low-x lines 
are excluded for the most metal-poor stars. 

To date, there is no general consensus about the mag- 
nitude of non-LTE corrections for iron, due to the in- 
completeness of the Fe model atom and the uncertainty 
about the rate of collision with the hydrogen atoms. 
Adopting several recipes for the calibration of the Sh 
parameter, the scaling-factor to correct the H I collision 
rate provided by the Steeenbock fc Holweger generalisa- 
tion dSteenbock fc Holwegerl 1 19841 ) of the Drawin formula 
(|Drawinlll968l . ll969l ). different authors provid e different non 
LTE corrections for iro n (see for instance iGratton et al.l 
ll999l;lGehr en et al.ll2004l . providing very diferent results). 

iMashonkina et all (|201ll ) analyzed a metal-poor giant 
star (slightly colder than those discussed here) by consid- 
ering the LTE case and the non LTE case with different 
efficiencies for the collision rate, and found that in case of 
Sh=0.1 (corresponding to a low efficiency of collisions with 



Li abundance in lower RGB stars 7 




5500 5000 4500 

T eH (K) 

Figure 2. Position of the target stars (grey points) in the 
T e ff— log g diagram, compared with the a-enhanced theoretical 
isochrones from the BaSTI database, for a reference age of 12.5 
Gyr and different metallicities (namely Z=0.0001, 0.0003, 0.0006 
and 0.001). The isochrones are plotted only up to the tip of the 
RGB, for sake of clarity. The location of the RGB bump is marked 
as a grey shaded region. The fainter shaded region marks the 
boundary beyond which A(Li) stops decreasing, i.e. attains its fi- 
nal value along the RGB . We also display two lines corresponding 
to A(Li) increased by 0.01 dex and 0.10 dex, (short-dashed and 
long-dashed, respectively). 

H I) a satisfactory excitation equilibrium is reached only by 
decreasing the temperature of ~ 80 K . 

Recently iBarklem et al.l (|201lh have convincingly 
shown that the Drawin formula lacks the necessary phys- 
ical ingredients to properly capture the quantum mechani- 
cal processes involved in excitation and ionisation, through 
collisions with hydrogen atoms. In fact in the cases for 
which such quantum mechanical computations have been 
performed the results differed by several orders of magni- 
tude with respect to the predictions based on the Drawin 
formula. This casts some serious doubts on the results of 
any NLTE computation which depends sensitively on the 
use of the Drawin formula and a suitable scaling factor. 

In our case the good agreement between our photomet- 
ric and spectroscopic T e ff, regardless of the metallicity, sug- 
gests that no significant non-LTE effects on iron are at work, 
at least for the lines of high excitation. 



4 THEORETICAL ANALYSIS 

Our theoretical analysis is based on a reference grid of 
stellar evolution models calculated with and without 
including the effect of atomic diff usion. More details about 
the code and input physics are in IPietrinferni et all (|2006t ) 
and iMucciarelli et all ()201ll ) . We have considered the a- 
enhan ced ([a/Fe]~0.4) metal mixture of IPietrinferni et al.l 
(|2006l ) and pairs of He and metal mass fractions (Y,Z) 



equal to (0.245,0.00001), (0.245,0.0001), (0.245,0.0003), 
(0.248,0. 002), (0.251, 0.004) that correspond to 
[Fe/H] = -3. 62, -2.62, -2.14, -1.31, and -1.01, respec- 
tively. The adopted AY/AZ=1.4 is the same as in 
IPietrinferni et al]|2006l ). For each metallicity we calculated 
the evolution of one stellar mass whose age along the 
RGB is equal to ~12.5 Gyr, and followed its evolution 
from the pre-MS to luminosities beyond the RGB bump. 
The range of mass values is comprised between ^0.80 
and ~0.86Mq, increasing with increasing Z. The initial 
A(Li) for all metallicities was fixed at the BBN predicted 
value of 2.72 dex. We have also calculated models with an 
initial Li abundance increased or decreased by one order of 
magnitude and found - as expected - that the amount of 
depletion after the first dredge-up (A (Li)) is insensitive to 
the exact value of the initial A(Li). 

Figure [3] displays the evolution of the surface Li abun- 
dance of models with initial [Fe/H] = — 1.01. The evolution 
without diffusion displays a 0.03 dex pre-MS depletion, fol- 
lowed by the MS evolution with unchanged surface A(Li). 
Dilution due to the deepening of the convective envelope at 
the end of the MS starts at T e ff ~5900 K, and continues 
until the first dredge-up is completed, at T e ff ~5100 K, 
when the model is already evolving along the RGB. As 
for the evolution with fully efficient atomic diffusion, the 
surface Li decreases during the MS and reaches a mini- 
mum around the turn off. After this point A(Li) at first 
increases, due to the deepening convection that returns to 
the surface part of the Li diffused outside the envelope along 
the MS. After a local maximum at T eff ~5850 K, A(Li) 
starts to drop, for the convective e nvelope reaches regions 
where Li has been bu rned (see, e.g. IDelivannis et alj|l99d : 
ISalaris fc Weissll200ll ). The surface Li abundance at the end 
of the dredge-up is only 0.07 dex lower when diffusion is 
fully efficient. 

Table [^displays A (Li) (in dex) for the metallicity range 
spanned by our calculations, and an age of 12.5 Gyr. An im- 
portant point to notice is that the [Fe/H] values reported 
in the first column are the initial values. In models without 
diffusion [Fe/H] at the surface stays constant all along the 
MS and is negligibly affected by the increase of surface He 
(and associate decrease of H) after the first dredge-up. How- 
ever, when diffusion is fully efficient, it decreases during the 
MS, because of the sinking of Fe and increase of H in the 
envelope. It is only after the first dredge-up that [Fe/H] is re- 
stored to within ~0. 02-0. 03 dex of the initial value. Working 
with lower RGB stars removes also the ambiguity between 
initial and actual [Fe/H] when diffusion is efficient. 

Overall, the total Li depletion from the pre-MS to the 
RGB shows a mild dependence on the initial metallicity. 
Metal rich models display a more efficient depletion, due es- 
sentially to deeper convective envelopes alo ng the RGB. This 
is in agree ment with previous resul t s by IDelivannis et al.l 
(|l990t ) and ISackmann fc Boothrovdl (| 19991 ). The values of 
A (Li) cover a range of ~0.15 dex over the [Fe/H] interval 
spanned by our calculations. The variation of A(Li) is non- 
linear, with an increase of A (Li) /A [Fe/H] with increasing 
[Fe/H] . 

We have also calculated models without diffusion but in- 
cluding - starting after the end of the MS - overshooting be- 
low the convective envelope. The extension of the overshoot- 
ing region has been parametrized in terms of the pressure 
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Figure 3. Evolution of the surface A(Li) as a function of T e ff 
for 0.855Mq, a-enhanced, [Fc/H]=— 1.01 stellar models without 
(solid line) and with (dashed line) the inclusion of atomic dif- 
fusion. The pre-MS, turn off (TO), subgiant branch (SGB) and 
RGB regions are marked in the figure. 



Table 2. Li abundance depletion (A(Li), in dex) along the RGB 
for the labelled [Fe/H] values. We report the results for mod- 
els without diffusion, with the inclusion of diffusion, and without 
diffusion but accounting for overshooting below the envelope con- 
vection, respectively (see text for details). 



[Fe/H] 


A (Li) (no diff) 


A(Li) (diff) 


A(Li) (oversh) 


-3.62 


1.28 


1.35 


1.29 


-2.62 


1.30 


1.37 


1.31 


-2.14 


1.33 


1.40 


1.34 


-1.31 


1.40 


1.46 


1.41 


-1.01 


1.44 


1.51 


1.50 



scale height (H p ) at the Schwarzschild boundary. We assume 
instantaneous mixing and radiati ve temperature gradien t in 
the overshooting region (see also iDelivannis et al. I ll990l , for 
similar calculations but a much larger overshooting region) . 
The need for these additional calculations stems from the 
results of analyses of the RGB bump luminosity in Galac- 
tic globular clusters. The predicted bump luminosities from 
models calculated with the same code and input physics 
adopted here, are too bri ght compared to o bservations, by 
^0.2 mag on average (see lCassisi et al.ll201ll , and references 
therein) although there is a hint of an increase of the discrep- 
ancy towards the lowest metallicities. In these additional cal- 
culations we have considered an overshooting distance equal 
to 0.35 Hp for the metal poor models up to [Fe/H]=— 2.14. 
This overshooting length lowers the RGB bump magnitude 
by ~0.40 mag. For models with [Fe/H]=-1.31 and -1.01 we 
have considered an overshooting distance equal to 0.10 H p , 
that lowers the bump level by ~0.2 mag. The correspond- 
ing values of A(Li) along the lower RGB are reported in 



Table 3. Estimates of the cosmological Li abundance (A(Li) ) 
in our sample of lower RGB Halo stars. We have considered three 
sets of stellar evolution models with different assumptions about 
the element transport mechanisms, and three T e ff scales for the 
A(Li) determinations (see text for details). 



Models A(Li) (A99) A(Li) (G09) A(Li) (spect) 

Standard 2.28 2.39 2.30 

Diffusion 2.35 2.46 2.37 

Overshooting 2.29 2.40 2.31 



Table [2] The effect of overshooting from the convective en- 
velope along the RGB - with our calibration - is negligi- 
ble (A(Li) increases by only 0.01 dex) with the exclusion 
of the model with [Fe/H] = -1.01 dex for which the differ- 
ence is of 0.06 dex. We have verified - with calculations for 
[Fe/H] = -3.32, -1.31 and -1.01 - that the effects of dif- 
fusion and overshooting (at least for our choice of the over- 
shooting extension) on A(Li) are additive, when we compute 
models including both effects. 

In addition, we have explored the effect of varying age, 
initial He content and mixing length, by calculating selected 
models at [Fc/H] = — 2.14. We found that increasing the age 
by 2 Gyr above our reference value of 12.5 Gyr decreases 
A(Li) by 0.02 dex. A decrease of the mixing length by 0.2H P 
below our solar calibrated value (2.01H P ) leaves A(Li) un- 
changed, while an increase of the initial He mass fraction Y 
by 0.01 increases the depletion by 0.01 dex only. 

4.1 Estimate of the primordial Li abundance in 
Population II stars 

As we have discussed in the section on the observational 
data, the measured values of the 12 C/ 13 C ratio and the po- 
sition in the T e / /-log g diagram place our star sample in the 
lower RGB phase, after the completion of the first dredge-up 
and before the onset of the mixing episode at the RGB bump 
luminosity. We have determined their primordial A (Li) - 
that we denote as A(Li) , to compare with the A(Li) Q pre- 
dicted by the BBN nucleosynthesis - for each of the three 
adopted T e j j scales, making use of a Monte Carlo technique, 
as detailed below. 

For each temperature scale we have first fixed the set of 
theoretical models. From the chosen set, we have produced 
a large (500 objects) sample of lower RGB synthetic A(Li) 
abundances, by drawing random values (flat distribution, 
consistent with the distribution in the observed sample) of 
[Fe/H] within the range spanned by our models, and deter- 
mining the corresponding depletion A (Li) by linear interpo- 
lation among the values listed in Tabled 

The value of A(Li) has then been determined as fol- 
lows. To each synthetic A(Li) we have added a constant 
that corresponds to a trial value for the initial Li abun- 
dance of the sample. Each element of the resulting ([Fe/H], 
A(Li)) synthetic sample has been then perturbed by a ran- 
dom Gaussian error. These error values have been calcu- 
lated taking the mean < a t s > spectroscopic error from 
the observations, perturbed by a Gaussian random distri- 
bution with standard deviation equal to the la dispersion 
of the spectroscopic errors around the value < a a t s >. The 
initial Li abundance of the synthetic sample is then varied 
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Figure 4. Comparison between a synthetic sample of A(Li) val- 
ues from standard models and A(Li) =2.30, and the observed 
A(Li) obtained employing the spectroscopic T e ff scale (see text 
for details). 

until the mean A(Li) of the synthetic sample matches the 
observed value. The corresponding initial Li-abundance is 
considered to be the best estimate of A(Li) Q . Finally, we 
have performed a Kolmogorov-Smirnov test to compare the 
synthetic distribution calculated with the best estimate of 
A(Li) with the observed one, and verified that, for all our 
choices of models and T e / / scale, the probability P that the 
two samples are not drawn from the same distribution is 
well below the 95% threshold. A comparison between syn- 
thetic and observed samples for a particular choice of the 
T e ff scale and stellar models is displayed in Fig. 2] 

Our final estimates of A(Li) are reported in Table 
Models that include atomic diffusion without overshooting 
predict A(Li) larger by 0.07 dex with respect to the stan- 
dard models, while models including only overshooting lead 
to Li abundance larger by only 0.01 dex when compared with 
the standard case (at least for the metallicity range covered 
by our targets). Inclusion of both diffusion and overshooting 
would increase A(Li) of 0.08 dex when compared with the 
standard models. 



5 LI ABUNDANCES IN LOWER RGB 
GLOBULAR CLUSTER STARS 

It is interesting to investigate also the surface A(Li) in lower 
RGB stars belonging to Galactic globular clusters, in com- 
parison with the values obtained for field stars. A note of 
caution is however warranted. The estimate and interpreta- 
tion of the initial Li abundance of Galactic globular clus- 
ters in terms of Halo primordial Li should be considered 
with caution, because of the occurrence of self-enrichment 
processes during the formation/evolution of the clusters, 
that are potentially able to affect the measured A(Li). A 



point to recall is that, in light of the currently accepted 
self-enrichment scenario for the evolution of the globular 
clusters, these systems lose a relevant fraction of their first 
generation stars - whose chemical abundance pattern is ex- 
pected to be consistent with the field Halo population - 
during their early evolution. Spectroscopic s tudies (see e.g. 
ICarretta et"aHl20ld : [ivlartell fc Grebdl20ld ) show that the 
large majority of field Halo stars have in fact an abundance 
pattern consistent with those of first generation stars in glob- 
ulars, hence both kind of objects are compatible in terms of 
A(Li). Only a very marginal fraction (less than 2%) of field 
halo stars is compatible with second generation cluster stars 
- that show the well known CN and ONa anticorrelation and 
have Li abundances potentially modified by nuclear process- 
ing. 

Lithium abundances are available (from the analysis 
of dwarfs and/or giants) for a few of the closest clus- 
ters. To date NGC 6752 (|Shen et all |2010| ) and 47 Tuc 
(|D'Orazi et al.l l20ld ) exhibit clear signatures of intrinsic 
star-to-star dispersion in A(Li), whereas M4 show a very 
high d egree of homogeneity in A(Li) at a given evolutionary 
phase (|Mucciarelli et al]|201ll ). NGC 6397 shows also sim- 
ilar homogeneous A(Li) abundances, with the detection of 
only 3 (out of 100) dwarf stars with A(Li)<2 (|Lind et al.l 
I2009D . 

The interpretation of these results is beyond the pur- 
pose of this work, and we discuss here only the Li abun- 
dance in giants stars of these clusters (namely NGC 6397, 
NGC 6752 and M4) for which high-resolution spectra of gi- 
ant stars are available, to provide a comparison with the Li 
abundances inferred in field RGB stars. 



5.1 NGC 6397 

We have analyzed 45 RGB stars (below the RGB bump 
level) in NGC 6397, retrieved from t he dataset of GI - 
RAFFE/FLAMES spectra discussed by iLind et all (|2009l ). 
All the observations are performed with the HR13 and HR15 
gratings. Due to the limited number of available iron lines, 
a full spectroscopic analysis cannot be performed and we 
discuss only the Li abundances inferred by employing the 
photometric T e // scales by A99 and G09. We employed 
the (J — K)o color, obtained from 2MASS photometry, cor- 
rected for reddening using the value of E(B-V) provided by 
iFerraro et al.l jl999l ). To reduce the scatter in the derived 
A(Li) distribution due to the photometric errors, we pro- 
jected each star along the fiducial line of the RGB sequence 
in the K-(J-K) plane (note that this technique, obviously 
applicable to stars members of a stellar cluster and not to 
field stars, is correct under the assumption that the broad- 
ening of the RGB is due to photometric er rors only). 

In their orig inal work lLind et al.l (120091 ) derive T e / / val- 
ues by_jKin£_thesyirUietic Stromgren (b-y) colors computed 
by lOnehag et all (|2009l ) employing MARCS model atmo- 
spheres, and find an average A(Li)= 1.13 dex (a= 0.09 
dex) for th e lowe r RGB sample. We find that the G09 and 
ILind et all l|2009h T e ff scales well match, with an average 
difference T°f f -T^ f y d09 = -13 K (cr= 15 K), while the differ- 
ence with the temperatures by A99 is Tffi - T e L /} ld09 = -117 
K (a= 14 K). The iron content is [Fe/H]= -2.08 dex (<r=0.09 
dex) with the G09 scale and -2.12 dex (<r=0.08 dex) with 



10 Mucciarelli et al. 



the A99 scale. We find an average A(Li)=1.00 dex (cr= 0.09 
dex) with the A99 scale and of A(Li)=1.09 dex (<r= 0.10 
dex) with the G09 scale. 

By interpolating in [Fe/H] among the values of Tab. 2, 
we derive an estimate for A(Li) of 2.33-2.42 dex with the 
model without diffusion (A99 and GH09 scales respectively) 
and 2.40-2.49 dex with the model with diffusion. 



5.2 NGC 6752 

Twenty-one RGB stars were observed with UVES within 
the ESO Large Program 65.L-0165 (PI: Grundahl). Con- 
cerning the Li abundance, only the qualitative behaviour 
of the obs erved EWs as a func tion of V magnitude was 
discussed (|Grundahl et al.l l2002h , without the explicit de- 
termination of the Li abundance. Here we consider the 12 
stars fainter than the RGB bump, for which the Li line is 
clearly detectable. We analyzed the spectra obtained with 
the CD#3 cross-disperser, deriving T e // both spectroscop- 
ically and photometrically, adopting the A99 and G09 cali- 
brations for the (J — K)o colour. The J and K s magnitudes 
are from the 2MASS da tabase, corrected for reddening us- 
ing the E(B-V) value bv lFerraro et all (|l999T ). Also for this 
cluster, we minimized the scatter in the final A(Li) values by 
employing (J — K)o colours obtained projecting the position 
of e ach star along the RGB fiducial line. 

iGrundahl et all 1)20021 ) derived the T e / / for their targets 
by means of the A99 calibration for the Stromgren (b-y) in- 
dex. We find a reasonable consistency with our spectroscopic 
and photometric T e // on the A99 scale (as for the case of 
field stars) with average differences T^fp - T °J ndahm2 = - 
35 K (<r= 45 K) and T$ - T°f f undahl02 = -43 K (cr= 25 
K), respectively. The offset with the T e ff by G09 is of 
T ?S B S - T°f f undahm2 = +65 K (er= 26 K). 

The average iron content is [Fe/H]= -1.68 dex (cr=0.07 
dex), when the spectroscopic temperatures are used and dif- 
ferent by only a few hundredths dex when the photometric 
T e ff are adopted. The average A(Li) is 0.83 dex [a= 0.15 
dex) with the spectroscopic temperatures (almost the same 
abundance is obtained with the A99 scale, while the G09 
scale provides A(Li)= 0.93 dex, a— 0.15 dex). Three stars 
display A(Li)~0.5-0.6 dex, lower by ~0.3-0.4 dex compared 
to the other stars in the sample. The remaining 9 stars pro- 
vide an average Li abundance A(Li)= 0.91 dex (<r= 0.07 
dex). Th e presence of these thr ee Li-poor stars (already iden- 
tified bv lGrundahl et alj|20b"3 ) is probably ascribable to the 
lithium variatio n in the cluster, as testified by the Li-Na 
anticorrelation (jPasquini et al.l [20051 ) and Li-0 correlation 
|Shen et al.ll2010r ) . These stars could belong to the second 
generation. 

By employing the A(Li) values listed in Tab. 2, the 
models without and with the inclusion of the diffusion pro- 
vide A(Li) estimates of 2.29-2.35, 2.18-2.24 and 2.19-2.25 
dex, by adopting, respectively, the T e // scale from GH09, 
A99 and from the excitation equilibrium. When we exclude 
the 3 Li-poor stars, the derived A(Li) increase by 0.08 dex. 

5.3 M4 

We summarized brie fly the results discussed in 
I Mucciarelli et all (j201lT ) about the metal-rich cluster 



M4, from the analysis of a sample of GIRAFFE spectra. 
The differential reddening that affects the field of view 
of M4 makes uncertain the abundances derived directly 
by adopting the photometric T e ff, due to the residual 
of the differential reddening correction. The atmospheric 
parameters for the RGB stars in the sample were therefore 
derived by projecting the position of each star along the 
stellar isochrone the best fit the observed colour-magnitude- 
diagram. This T e ff scale is in nice agreement w it h the 
spectroscopic T e // values inferred by iMarino et all l|2008T ) 
for the stars in common. 

The derived Li abundance in lower RGB stars of M4 
is A(Li)=0.92 dex, that leads to A(Li)o=2.35 dex (models 
without diffusion) and 2.40 dex (models with diffusion). 



6 SUMMARY AND CONCLUSIONS 

We have discussed the use of Li abundances measured in 
Population II lower RGB stars as an independent, reliable 
and robust diagnostic of the initial Li abundance in the 
Galactic Halo. Surface abundances in giant stars fainter than 
the RGB bump are sensitive to the total Li content left at 
the end of the MS phase, and are very weakly affected by 
atomic diffusion during MS. Also, the predicted A(Li) in 
these objects is basically insensitive to the mixing length cal- 
ibration, the precise stellar ages, initial He abundances, and 
realistic estimates of the overshooting extension below the 
Schwarzschild boundary of the convective envelope. Chemi- 
cal abundance measurements in giant stars also suggest that 
any additional element transport along the RGB is very 
likely inefficient in this phase. Overall, our analysis reveals 
that the predicted Li depletion A (Li) along the lower RGB 
is robust in terms of theoretical interpretation. 

The values of A(Li)o inferred from our sample of lower 
RGB Halo field stars range from 2.28 (obtained with the 
A99 scale and without the inclusion of atomic diffusion) to 
2.46 (when the G09 T e ff scale is used, together with models 
including atomic diffusion). Inclusion of overshooting from 
the RGB convective boundary would increase both limits by 
only 0.01 dex, while variations of age, initial He abundance 
and mixing length parameter provide changes of a few hun- 
dredths of dex or less in terms of A (Li). When a T e // scale 
is adopted, the effect of fully efficient atomic diffusion on 
the A(Li)o estimate is by at most 0.07 dex. 

The discrepancy with A(Li)o predicted by BBN calcu- 
lations thus remains, reconfirmed by the robustness of our 
estimate. The analysis performed on the lower RGB stars in 
three Galactic globular clusters confirms similar values for 
A(Li)o, although in general the possible occurrence of self- 
enrichment processes in these objects has to be considered. 

There are several discussions in the literature about 
how to solve this discrepancy with BBN results. Ideas in- 
volve a first generation of stars that has processed and ef- 
ficiently depleted lithium in a substantial frac tion of the 
early Halo baryonic matter (|Piau et al.l l2006f ). modifica- 
tions to t he BBN considering the decay of unstable par- 
ticles (e.g. ICvburt et al.l 120101 ). modifications to the reac- 
tion cross sections for th e Li production during BBN (e.g. 
IChakrabortv et "aT1l201lf ). or 'deep' turbulent mixing during 
the MS that connect s the convective envelope to inner, Li- 
depleted regions (e.g. iRichard et al.|[2005l ; iMucciarelli et all 
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I2011I ). Our result for RGB stars provides an additional ro- 
bust constraint on the efficiency of this turbulent mixing. 
Independently of its parametrization, during the MS any 
additional element transport needs to bring into the Li 
burning region an amount of initial lithium A;, urn (Li)=0.3- 
0.4 dex, in order to eliminate the discrepancy with BBN cal- 
cul ations (the turbulenc e model denoted as T6.25 or T6.28 
bv lRichard etU] (120051) appear to burn approximately the 
right amount of Li in models with [Fe/H]= -2.31). The value 
of A;,„ rn (Li) appears to be roughly the same for both Halo 
field and globular cluster stars. 

The use of lower RGB stars will allow estimates of the 
Li content in stellar populations more distant that those 
usually observed to investigate the Spite Plateau. The 
obvious benefits are: (i) an enlarged sample of field stars 
and clusters to study the primordial Li abundance within 
the Galaxy, and (ii) a better chance to assess whether a 
"cosmological lithium problem" exists also in extragalactic 
systems with a different origin and a different star formation 
history (see iMonaco et al.l l201ll . for a first study of the 
initial Li in u Cen, widely considered as the remnant of a 
dwarf galaxy accreted by the Milky Way ). For instance, 
the current generation of high resolution spectrographs 
mounted on 8 metre-class telescopes allows to reach down 
to one magnitude fainter than the RGB bump level for stars 
in M 54 and the old stellar population of Sagittarius dwarf 
galaxy. In principle, also stars at the RGB bump level in 
old clusters belonging to the Large Magellanic Cloud can 
be reached, by using, for example, the VLT spectrograph 
X-SHOOTER; this possibility (even if expensive in terms of 
observing time due to the faintness of the targets) will open 
a new perspective in the investigation of the Li problem. 
Finally, with the advent of 30 metre-class telescopes (as the 
E-ELT) spectroscopy of lower RGB stars will provide new 
constraints to the primordial Li abundance in even more 
distant systems. 
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